Science Advances

Bone remodeling is a tightly regulated physiological process that maintains bone health through coordinated action of bone-resorbing osteoclasts and bone-forming osteoblasts. Disruption of this balance, such as the one induced by estrogen decline after menopause, results in bone loss and osteoporosis. Genetic factors play an important role in determining bone mineral density (BMD) loss over time. However, translating genetic associations into individualized risk prediction remains challenging due to small effect size of individuals variants and non-linear interactions within the bone remodeling unit. Here, we present a bone cell population dynamics model that includes major regulatory pathways, such as the RANK/RANKL/OPG axis, Wnt signaling, and hormonal regulation by estrogen, parathyroid hormone, and TGF-{beta}. We calibrate the model on clinical data from healthy postmenopausal women, and women with reduced BMD undergoing anti-osteoporotic therapy. The calibrated model captures healthy BMD decline in postmenopausal women and therapeutic response to anti-osteoporotic medications. We mechanistically incorporate the effect of 22 variants across 8 genes involved in bone remodeling and simulate BMD trajectories in 1,000 virtual subjects differing by ancestry and genetic makeup. The median predicted 5-year BMD loss was 3.57% (95% prediction interval: 1.31-5.24), consistent with the values reported in the literature. The virtual individuals with African ancestry were predicted to experience the highest average 5-year BMD loss. The strongest genetic risk factors for bone loss were predicted to be CYP19A1 rs727479 and OPG rs3102735, while LRP5 rs11228240 emerged as a protective factor that could partially counteract the detrimental effects of other variants. Several epistatic effects were observed in the genetic interaction analysis. Mechanistically, our model suggested that estrogen exerts its effect on bone remodeling primarily by modulating osteoclast apoptosis. Overall, this framework demonstrates a proof-of-concept for integration of genetic risk factors into mechanistic models of disease and can be extended to other conditions with polygenic inheritance.

The recent availability of at-home menstrual cycle tracking technology has created opportunities for personalized assessment of reproductive health, alongside improved characterization of hormone patterns in women with and without reproductive disorders such as polyendocrine metabolic ovarian syndrome (PMOS), which affects approximately 10% of reproductive-age women. In this study, we leverage self-tracked urinary hormone data to develop an autoregressive Hidden Markov model (arHMM) that maps cycle days to physiologically meaningful phases based on hormone trajectories. By modeling day-to-day hormonal dynamics rather than absolute hormone levels, and allowing variable phase durations, this approach accommodates substantial variability in menstrual cycles, thereby enabling meaningful comparisons within and between individuals. Across more than 3800 cycles from over 1100 individuals, we find that arHMM-derived phases reproduce expected hormonal patterns within follicular, periovulatory, and luteal phases, and that phase-based timing for hormone testing outperforms conventional cycle day-based testing in capturing the luteinizing hormone surge and post-ovulatory progesterone rise, highlighting limitations of fixed-day clinical protocols. We identify phase-specific differences between healthy controls and individuals with self-reported PMOS, including lower luteinizing hormone in the periovulatory phase, and reduced luteal-phase progesterone levels in PMOS. Furthermore, features derived from arHMM phase assignments enable classification of PMOS status with ~78% accuracy, demonstrating the potential of this approach for non-invasive PMOS screening.

Accurate assessment of body composition is important to risk stratification and management of metabolic, musculoskeletal, and aging-related diseases, yet reference modalities such as Dual-energy X-ray absorptiometry (DXA) are costly and impractical for frequent monitoring. Commodity 3D body scans offer a low-cost, radiation-free alternative, but extracting meaningful and predictive shape features from scans remains challenging due to nonuniform point density, variable body size and cross-device differences. We introduce BodyMAE, a self-supervised, surface-area aware masked autoencoder for metric-scale 3D body scans. The pipeline integrates area-adjusted sampling, a long-range focused encoder, and a lightweight decoder regularized to promote locally uniform reconstructions. Trained and evaluated on 917 paired 3D body scans paired with clinical DXA reports, BodyMAE achieves strong accuracy on fat percentage (root-mean-square error (RMSE) 3.825 percentage points, R^2 0.908), fat mass (RMSE 3.694 kg, R^2 0.968), and lean mass (RMSE 3.608 kg, R^2 0.901), with competitive performance on bone mineral content (RMSE 0.284 kg, R^2 0.754).We also assess feature stability across pretrained baselines, finding higher retrieval accuracy for our representations (Top-1 90.131%). These results indicate that combining metric-aware sampling, long-range relational encoding, and local geometric regularization enables accurate body composition estimation from 3D body scans, as validated by comparisons to DXA-derived measurements.

Mood and anxiety disorders emerge predominantly in adolescence, yet they are usually identified only once symptoms have consolidated, when intervention can only be reactive. A marker that registers the loss of healthy brain function before symptoms crystallise would allow earlier and more targeted treatment, much as caged canaries once warned miners of danger before it became apparent. Here we report such a marker using a single baseline resting-state functional MRI scan in 150 adolescents in the Human Connectome Project Boston Adolescent Neuroimaging of Depression and Anxiety (HCP BANDA) cohort, allowing us to prospectively predict depression and anxiety symptoms one year later in held-out participants at r = 0.60, substantially above the effect-size ceiling reported for functional connectivity in the same data. The marker is not computed from raw functional connectivity but read out from a whole-brain generative model fitted to each individual's dynamics, which gives access to interference structure that covariance-based features cannot represent. The regions driving the prediction, including precuneus, ventromedial prefrontal and anterior cingulate cortices, are among those previously implicated in internalising disorders, and the same signature tracks cognitive variation in healthy participants and is mechanistically linked to the efficiency of task-related computation. These findings establish a mechanistically interpretable and prospectively predictive marker of adolescent mental health and define a clear path towards external validation and clinical use.

Postoperative resilience varies widely among older adults, yet the biological drivers of recovery remain unclear. We evaluated whether preoperative immune profiles, measured in plasma and through ex vivo whole blood stimulation, predict resilience to the acute stress of total knee arthroplasty. A total of 152 adults (greater or equal to 60 years) in the PRIME KNEE cohort underwent elective total knee arthroplasty and had available blood samples for measurement of 45 immune biomarkers, quantified in plasma and in whole blood stimulated ex vivo for 24 hours with lipopolysaccharide (LPS) or influenza antigen (FLU). Resilience was assessed using Expected Recovery Differential (ERD) and Resilience Trajectory (RT) across pain severity, pain interference, lower extremity physical activities of daily living (LE PADLs), and step counts. An exploratory stability selection framework using LASSO identified biomarker predictors of postoperative outcomes. Plasma and stimulated biomarkers showed broadly similar predictive performance. A shared set of biomarkers, including LBP, leptin, TNFR1, CD30, and LIF, was consistently selected across models. Immune predictors explained ~12-24% of the variance in resilience outcomes. Distinct immune signatures emerged for pain versus functional recovery: pain related predictors mapped to local inflammatory and neuroimmune pathways, whereas function related predictors reflected systemic inflammatory load and cytokine signaling. Preoperative immune biomarkers, whether measured in plasma or after ex vivo stimulation, capture meaningful variance in postoperative resilience. The divergence between pain related and function related immune signatures highlights biologically distinct pathways underlying different dimensions of recovery and supports further development of immune based perioperative risk assessment.

Biological systems evolve as continuous dynamical processes, but at organ-scale and across human lifespans they are rarely observed longitudinally--population data typically exist instead as sparse, cross-sectional snapshots. Inferring lifespan dynamics from such data requires methods distinct from those used at cellular and tissue scales where dense observations are accessible. We address this problem in the thoracic aorta, where surgical decisions currently rest on static, age- and sex-agnostic diameter thresholds that reduce three-dimensional morphology to a single scalar. Treating normal aortic morphology as a stochastic dynamical system, we pose a continuous-time drift-diffusion process in a two-coordinate state space of normalized surface area (A) and normalized fluctuation in integrated Gaussian curvature ({delta} K), and fit closed-form solutions of the Fokker-Planck equation by maximum likelihood to a sex-balanced, age-uniform cohort spanning infancy to age 99. Inter-individual variability is treated as a fitted diffusion parameter rather than as residual scatter, which is distinct from prior normative studies that report variability as scatter around a regression line. The framework identifies two growth regimes for aortic size (childhood expansion followed by persistent adult growth, with adult males growing approximately 70% faster than adult females) and a single dynamical regime for aortic shape, with heteroscedastic variability accumulating at a rate comparable to the mean drift over the lifespan. Applied to independent cohorts of acute and chronic thoracic aortic dissections, the multivariate model identifies over 95% as statistical outliers via Mahalanobis distance, consistently outperforming either coordinate alone. The same probabilistic envelope that describes normal aging thus defines a baseline against which disease can be detected, supporting a shift toward dynamic, age- and sex-aware assessment of thoracic aortic pathology.

Introduction: Acute low back pain is a leading cause of disability worldwide. Clinical guidelines recommend non-pharmacological therapies as first-line treatment and advise caution with opioid prescribing. However pharmacological therapies, including opioids and gabapentinoids, remain commonly used. The comparative risks of subsequent opioid use disorder (OUD) and overdose diagnosis associated with initial treatment modality in large, real-world populations is not well characterized. We estimated the incidence of new-onset OUD and overdose diagnosis among opioid-naive, Medicaid-insured adults with newly diagnosed acute low back pain and estimated the association between initial treatment modalities and subsequent OUD and overdose diagnosis risk. Methods: We conducted a retrospective cohort study using Medicaid T-MSIS Analytic files from 25 states (2016-2019). We identified opioid-naive adults with a new diagnosis of acute low back pain who initiated pharmacologic or non-pharmacologic treatment within 1 month of diagnosis. The primary outcome was incident OUD and overdose diagnosis (based on diagnosis codes in claims) during follow-up. Associations between initial treatment modality and OUD and overdose diagnosis risk were estimated using a non-parametric, doubly robust estimator to adjust for measured confounding. Results: The cohort included 525,002 opioid-naive adults initiating treatment for low back pain. The cumulative incidence of OUD and overdose diagnosis was 1.5% and 2.4% at 7 and 13 months, respectively. Compared to non-use, use of gabapentinoids during the first month of treatment was associated with the highest relative risk (increasing risk) by 130.1%, 95% confidence interval (CI): 117.8%, 142.3%), the second-highest relative risk was estimated for higher-dose opioids, defined as > 50 daily Morphine Milligram Equivalents (MME) (118.1%, 95% CI: 99.2%, 137.0%). Lower-dose, short-duration opioids ([≤] 50 MME, [≤] 7 days) were also associated with elevated risk, though substantially smaller in magnitude (20.8%, 95% CI: 13.8%, 27.9%). In contrast, non-pharmacologic, non-interventional therapies were associated with reduced OUD and overdose diagnosis risk, with physical therapy demonstrating the largest relative reduction of 34.0% (95% CI: -40.9%, -27.1%). Discussion: In opioid-naive Medicaid patients with acute low back pain, initial non-pharmacologic treatment was associated with reduced OUD and overdose diagnosis risk. Gabapentinoids and opioids were each associated with increased risk; for opioids, the degree of risk increased with higher doses and durations. These results support guideline recommendations favoring non-pharmacologic treatment as first-line therapy and indicate the importance of cautious prescribing when pharmacologic treatment is considered.

Metabolic dysfunction-associated steatotic liver disease (MASLD) is rapidly increasing worldwide, yet effective targeted therapies remain limited. To better understand the molecular mechanisms underlying MASLD, we performed an integrated proteogenomic analysis of human liver tissue. Using mass spectrometry, we quantified 2,744 proteins in 504 liver biopsies from the Quebec Obesity Biobank and examined changes across disease stages. To investigate causality, we integrated liver proteomics with RNA sequencing and genome-wide genotyping to map thousands of protein quantitative trait loci (pQTLs) and expression quantitative trait loci (eQTLs). These molecular data were combined with summary statistics from a meta-analysis of genome-wide association studies including 16,532 MASLD cases and 1,240,188 controls. Mendelian randomization and genetic colocalization analyses revealed that most proteins differentially expressed across MASLD stages were not causally implicated in disease risk, whereas several genetically predicted liver proteins showed evidence of causal effects. Among these, higher hepatic levels of the MTARC1 protein were causally associated with MASLD and hepatic fat accumulation. Phenome-wide analyses suggested that MTARC1 inhibition may reduce the risk of cirrhosis, hepatocellular carcinoma, and cholelithiasis while improving lipid profiles. Notably, the causal MTARC1 variant influenced liver protein levels but not gene expression. Genetic analyses also identified ERLIN1 and HSD17B13 as potential therapeutic targets. In contrast, eQTLs and pQTLs at other loci such as GCKR showed opposite effects on MASLD risk. These findings highlight the importance of integrating tissue proteomics with human genetics to distinguish biomarkers from causal drivers and to identify promising therapeutic targets for MASLD.

Acute respiratory distress syndrome (ARDS) is a devastating complication of respiratory infections; however, the biological mechanisms that initiate its onset are poorly defined. Here we show that TNFRSF13B polymorphisms increase the risk of ARDS following SARS-CoV-2 infection up to 7.4-fold compared to the WT genotype. The increased risk was not due to immune-deficiency or impaired virus neutralization. On the contrary, TNFRSF13B mutant subjects mounted better antibody neutralization compared to subjects with WT TNFRSF13B. However, IgG from subjects expressing TNFRSF13B variants had less sialic acid, terminal galactose, and fucose than IgG from subjects with a WT genotype. Moreover, IgG from TNFRSF13B mutant subjects exhibited increased recruitment of complement factors. Thus, besides well-known actions governing plasma cell differentiation, TNFRSF13B impacts both affinity maturation and effector functions of IgG in ways that independently govern complement activation controlling inflammatory responses known to trigger ARDS.

Metabolic diseases such as type 2 diabetes (T2D) arise through complex interactions between physiological, molecular, and environmental processes. Clinical traits including age, sex, adiposity, and glycaemic status are strongly associated with disease risk and progression, yet most molecular studies examine these factors independently and assume relatively static molecular regulation. Consequently, how physiological state dynamically reshapes molecular organisation across omics layers remains poorly understood. Here, we integrated transcriptomic, proteomic, metabolomic, and genetic data from 3,027 individuals in the IMI DIRECT cohort to characterise the joint molecular effects of age, sex, body mass index (BMI), and glycated haemoglobin (HbA1c). We identified widespread associations between these traits and molecular phenotypes. However, interaction analyses revealed a more complex context-dependent regulation, showing that the molecular effect of one trait frequently depends on the state of another, with sex-specific effects of age being more prominent. We also investigated relationships between different types of molecular phenotypes and how these relationships are modulated by metabolic disease relevant traits, demonstrating that cross-omic molecular coordination is itself dynamically remodelled by physiological and metabolic state. Probabilistic causal inference identified a directionally structured network of age-associated molecules, revealing pathways through which age effects propagate across omics layers, showcased in the example of the mTOR signalling pathway. Integration of this directed network with genetic colocalisation analyses also identified a sub-network relevant for T2D. Collectively, our findings demonstrate that metabolic disease relevant traits not only independently influence molecular phenotype abundance but also jointly reshape the directional organisation of cross-omic molecular networks. These results support a model in which metabolic disease susceptibility emerges through dynamic rewiring of interconnected molecular systems and provide a framework for context-dependent biomarker discovery, disease stratification, and precision metabolic medicine.

Educational attainment-related polygenic scores have been implicated in autism spectrum disorder (ASD), but how parental polygenic scores shape offspring phenotypes remains unclear. Using genotyping and exome-sequencing data from 142,357 individuals (55,252 ASD cases) in a large ASD cohort, we dissected the direct and indirect genetic effects of educational attainment-related polygenic scores on ASD phenotypes. Trio-model analyses showed that parental polygenic scores for educational attainment (PGSEA ) were associated with milder core ASD symptoms, including social deficits and repetitive behaviors, predominantly through indirect genetic effects, whereas their associations with comorbidities were driven predominantly by direct genetic effects. PGSEA was also significantly negatively associated with rare variant burden and prenatal factors, although these factors contributed largely independently to most phenotypes. Adjustment for full-scale intelligence quotient (FSIQ) and socioeconomic status (SES) partially attenuated the indirect effects of PGSEA on offspring phenotypes. Finally, higher parental PGSEA was associated with later age at diagnosis in offspring, partly through its protective effects on ASD phenotypes. These findings indicate that indirect genetic effects of parentalPGSEA contribute substantially to phenotypic variation in ASD and highlight family-mediated pathways as an important component of ASD heterogeneity.

Background. The emergency department in the United States of America functions as a residual access point for healthcare and social services for populations including rural communities, the uninsured, mental health and addiction patients, and the unhoused. The workforce variable that determines unit function (experience density, the concentration of accumulated clinical judgment within a unit workforce) is not measured in hospital accounting systems. Objective. To document workforce composition changes in U.S. emergency nursing across the 2018 and 2022 cycles of the National Sample Survey of Registered Nurses (NSSRN), and to specify falsifiable predictions for the 2026 cycle. Methods. We analyzed NSSRN public-use files using a four-way ED definition extending Castner et al. (2024) and a hospital-bedside-restricted comparator. Variance estimation used jackknife replicate weights for 2018 and Successive Differences Replication for 2022. Burnout was operationalized using the Norful et al. (2023) leaving-reasons proxy across cycles, with sensitivity analysis using the 2022 direct burnout item. Results. A 15-year trajectory (2008-2022) documents progressive experience-density compression: the ED's 15+ year veteran cohort fell from 41.9% to 28.0% over the decade preceding the pandemic, a loss of nearly a third of the senior cohort and a 19.6% decline in mean experience density, before recovering modestly to 33.3% as veteran nurses remained through the pandemic acute phase, leaving the ED as the youngest hospital setting throughout. Hospital non-ED bedside nurses lost senior tenure between cycles (mean 15.65[-&gt;]14.06 years since first licensure; 15+ year share 43.5%[-&gt;]38.7%), while ED nurses retained their senior tail (mean 11.60[-&gt;]12.58). Burnout endorsement rose sharply in both populations (non-ED 27.3%[-&gt;]46.0%; ED 34.2%[-&gt;]61.2%), with the ED-vs-non-ED gap more than doubling. Controlling for tenure, ED status was not independently associated with burnout in 2018 (OR 1.15, 95% CI 0.83-1.59) but was strongly associated in 2022 (OR 1.92, 95% CI 1.44-2.55; p<.001). The direct burnout item showed a parallel pattern (OR 2.92, 95% CI 1.62-5.28). Conclusions. A pandemic-era setting-specific burnout effect emerged in emergency nursing that workforce-composition controls cannot explain. The 2022 cycle establishes a pre-exit baseline against which the 2026 NSSRN will serve as the falsifiable test of post-Omicron veteran exit. Nursing pipeline replacement lag exceeds the interval before 2026 data arrives; the consequences of inaction fall on populations dependent on ED-based residual access.

Pathogenic variants in ATP13A2, which encodes an endolysosomal polyamine exporter, cause Kufor-Rakeb syndrome and are associated with early-onset parkinsonism and related neurodegenerative disorders, however, the mechanisms by which ATP13A2 dysfunction drives disease remain incompletely defined. In Atp13a2 knockout mice, we identified an early, transient reduction in brain polyamines that precedes overt gliosis and behavioural abnormalities. Pharmacological polyamine depletion exacerbates phenotypes, whereas oral supplementation of spermidine, but not spermine, rescues parkinsonian symptoms establishing metabolic polyamine deficiency as a pathogenic driver. Mechanistically, spermidine counteracts microglia lysosomal dysfunction in the brain and exerts mitochondrial antioxidant and anti-inflammatory effects in primary mouse microglia, thereby improving neuronal integrity. In the absence of Atp13a2, microglial spermidine import relies on the related polyamine transporter Atp13a3. Importantly, these findings translate to human systems, whereby spermidine attenuates inflammation in ATP13A2-deficient human differentiated microglia, while postmortem ATP13A2-deficient brain analysis confirms increased microglia reactivity. Spermidine also rescues motor deficits and dopaminergic neuron loss in ATP13A2-deficient Drosophila and other fly parkinsonism models. Together, these findings identify early polyamine dysregulation as a mechanistic contributor to ATP13A2-associated parkinsonism and nominate spermidine supplementation as a potential therapeutic strategy for ATP13A2-driven pathology and possibly a broader range of parkinsonian sub-types.

Adapting ongoing gait patterns to environmental challenges is essential for safe navigation through the environment. Impairment of gait adaptation is common in many neurodegenerative disorders, such as Parkinson's disease (PD), where it hampers mobility and limits quality of life. The neural control of gait adaptation remains largely unclear, thereby limiting the development of targeted treatments, such as deep brain stimulation of the subthalamic nucleus (STN-DBS). We integrated clinical, kinematic, brain metabolic imaging, and electrophysiological data, obtained during a fully immersive virtual reality overground walking task, to characterize the neural underpinnings of gait adaptation performance during dynamic obstacle avoidance and its improvement with STN-DBS. Movement kinematics, brain oscillatory activity, and metabolic activation were simultaneously acquired in 12 patients with PD during rest and gait adaptation, under active or paused STN-DBS, using inertial measurement units, electroencephalography, and three separate [18F]fluorodeoxyglucose positron emission tomography scans. Eight age-matched healthy subjects completed the same task for comparative kinematic analyses. All patients showed significant clinical improvement with STN-DBS. During the gait adaptation task with paused stimulation, patients exhibited increased metabolic activity in the cerebellum and sensorimotor cortex. Active STN-DBS selectively enhanced thalamic and superior frontal gyrus (SFG) metabolism, while concomitantly reducing cerebellar uptake. Right-lateralized SFG metabolism correlated with gait adaptation performance, with DBS-driven shifts toward greater right SFG activity predicting the magnitude of gait adaptation improvement. This correlation was independent of baseline asymmetry in clinical impairment, electrode placement, or structural connectivity to the SFG. Of note, STN-DBS amplitude asymmetry emerged as an independent predictor of right-lateralization of SFG metabolism. EEG recordings confirmed this lateralized network modulation, with theta-band asymmetry paralleling PET findings. Our findings identify a lateralized thalamo-cortical network supporting gait adaptation in PD and highlight a distinctive role for the SFG. We further show that effective STN-DBS acts as a lateralized regulator, dynamically rebalancing cortico-thalamic circuits to support context-appropriate gait control. The observed right-hemispheric lateralization may foster novel image-guided programming strategies to enhance the consistency and effectiveness of gait control in PD.

Chronic kidney disease is characterized by decreased glomerular filtration rate (eGFR, estimated from serum creatinine or cystatin C) or increased urinary albumin-to-creatinine-ratio (UACR). Genome-wide association studies provided the genetic make-up of these traits, but their overlap remained largely unknown. Our multi-trait GWAS (N=1M) identified 812 signals and multi-trait fine-mapping sharpened the identification of likely causal variants. Of 333 signals classified for filtration function or albuminuria, only 11 overlapped. Their effects on eGFR and UACR were directionally concordant, dominated by eGFR and independent of HbA1c or mean arterial pressure. Mapped genes pinpointed mechanisms related to glomerular filtration area (SHROOM3, EPB41L5) and sodium-mediated intraglomerular pressure (NRBP1, DPEP1/CHMP1A). Genetics of fluid intake resulted in shadow effects on UACR without albumin leakage into urine. Our multi-trait approach sharpened the identification of likely causal genes for kidney traits, demonstrated largely distinct genetics for filtration function versus albuminuria, and provided new biological insights into the overlap.

Gene therapy is rapidly emerging as a transformative treatment for monogenic neurological disorders, including pediatric movement disorders such as aromatic L-amino acid decarboxylase (AADC) deficiency. However, its success critically depends on defining target cells and windows for therapeutic intervention. Here, we present an open-access single-nucleus transcriptomic atlas of the human basal ganglia spanning a therapy-relevant window from second/third trimester to the perinatal period and adulthood. Across 35,755 nuclei, we identify major (non-)neuronal cell types, retrace developmental trajectories, and characterize gene-regulatory networks. We identify so far unrecognized human-specific expression of key neuronal signaling genes, including GNAO1 and ADCY5, and discuss the implications for targeted gene replacement therapies. Unexpectedly, we found that the Huntingtin gene (HTT) is already expressed during prenatal stages of human brain development, supporting a previously proposed neurodevelopmental component of Huntington's disease, which should be considered in diagnostic and therapeutic strategies. Moreover, FOXG1 expression and regulon activity are predominantly located in a prenatal time window, suggesting constraints on the effectiveness of postnatal interventions. Our findings highlight the importance of datasets capturing human brain development in real time and provide a publicly available resource to guide precision gene therapy strategies in the future.

Post-surgical adverse outcomes, including mortality, intensive care readmission, and complications, remain major challenges for clinical decision-making. Existing machine learning approaches focus on outcome prediction while operating as opaque systems, limiting clinical trust and the translation of predictions into treatment decisions, and many clinical studies rely on synthetic data in which shared intermediate variables create circular dependencies between inputs and targets that compromise reported performance. We aimed to develop an explainable multimodal architecture and a rigorous evaluation methodology that address these gaps. We designed a two-stage architecture integrating supervised deep learning for risk prediction with conservative Q-learning for action recommendation. The first stage uses five modality-specific encoders for structured records, physiological time-series, chest radiographs, clinical notes, and surgical metadata, unified through cross-modal attention into a shared patient-state representation. The second stage applies offline reinforcement learning to recommend clinical actions while preventing value overestimation. We formally characterized a target-leakage flaw in synthetic pipelines and propose a real-data methodology using a verified clinical database, with event-censored temporal separation and uncertainty-weighted per-task training. Component-level behavior was validated on a controlled synthetic benchmark, demonstrating that the architecture functions as designed without claiming clinical validity. The cross-modal attention and risk-prediction components behaved as expected, whereas the offline reinforcement learning stage did not converge on the benchmark, indicating that value estimation requires further investigation on real clinical data. The architecture provides dual-level explainability through attention visualization and value decomposition, contributing a deployable design, a formal methodological critique of synthetic-data practices, and a complete framework for clinically valid evaluation.

Stimulant use disorder (StimUD) is a significant public health problem, but genetic studies have been limited by small sample sizes. We conducted genome-wide association studies (GWAS) of StimUD in the Million Veteran Program (MVP) and All of Us (AOU), followed by meta-analysis with FinnGen and 10 additional datasets, for a total of 709,369 individuals (Ncases=33,977, Ncontrols=675,392) in four broad ancestry groups: European (EUR) (Ncases=22,564, Ncontrols=624,672), African (AFR) (Ncases=7,574, Ncontrols=34,189), Admixed American (AMR) (Ncases=3,657, Ncontrols=15,698), and East Asian (EAS) (Ncases=182, Ncontrols=833). Population-specific SNP heritability was 6.1% in EUR and 2.4% in AFR. We discovered a total of 19 genome-wide-significant loci, six in EUR, including DRD2*rs5794864, P=7.32E-10, one in AFR, five in a multi-ancestry meta-analysis, including CHRNA5*rs55781567, P=3.27E-9, two in a male-only meta-analysis, including FTO*rs8057044, P=9.50E10-9, and five in a meta-analysis of sex-stratified results. In a hold-out AOU subsample (NEUR=18,841, NAFR=12,263, NAMR=9,739), ancestry-specific polygenic risk scores were significantly associated with StimUD in EUR (OR=3.28, 95% confidence interval (CI)=2.89-3.71) and AMR (OR=2.01, 95% CI=1.71-2.37). Transcriptome-wide association studies, fine-mapping, and colocalization analyses prioritized additional genes (e.g., GPX1, BSN). Genetic correlation, Mendelian randomization, and causal mixture analyses revealed relationships with other substance use and use disorder phenotypes, including cannabis use disorder (rg=0.94, P=5.43E-237) and opioid use disorder (rg=1.01, P=4.40E-107), and other psychiatric traits, including anxiety, depression, neuroticism, and attention-deficit/hyperactivity disorder. This is the first well-powered GWAS of StimUD, and it offers significant insights into disease biology.

Severe preeclampsia (sPE) is a major cause of maternal and fetal morbidity worldwide, yet its placental molecular heterogeneity remains poorly defined by current clinical diagnosis. To resolve the molecular architecture of sPE, here we integrated DNA methylation and proteomic profiling from a multi-ethnical cohort of 444 placentas from the Hawaiian Biorepository (HiBR), including 169 sPE cases, matched preterm controls and full-term controls. To address cellular heterogeneity in bulk placental tissue, we developed HOMED (Hierarchically Optimized Methylation Deconvolution), a single-cell-guided hierarchical framework for inferring placental cell-type composition from DNA methylation data. HOMED-adjusted integrative analyses identified extensive subtype-specific alterations involving hypoxia, angiogenesis, immune activation, trophoblast differentiation and metabolic remodeling. Molecular stratification revealed two reproducible sPE subtypes with divergent placental aging trajectories. One subtype exhibited a pre-mature placental state marked by accelerated placental aging, whereas the other displayed slower accelerated placental aging but a substantially increased risk of small-for-gestational-age birth (P = 0.028). These subtypes were independently replicated across six external cohorts and further supported by proteomic signatures achieving a classification accuracy of 0.88. Integrative epigenomic and proteomic analyses linked the growth-restricted subtype to hypoxia-associated glycolytic remodeling, suggesting distinct pathogenic mechanisms underlying clinically diagnosed sPE. Together, our findings redefine severe preeclampsia as a biologically heterogeneous placental disorder composed of molecularly distinct subtypes with divergent aging trajectories and fetal growth outcomes, providing a framework for mechanism-based stratification and precision obstetric medicine.

Whole-genome sequencing comprehensively captures coding, non-coding and structural variation in families with suspected inherited disorders, yet its clinical utility remains constrained by an interpretation bottleneck: selecting a handful of relevant variants from millions of candidates. Current rule-based pipelines, anchored in ACMG/AMP criteria, excel at identifying highly penetrant Mendelian alleles but frequently miss variants of low-to-moderate penetrance, non-coding alterations and germline-somatic interactions. Here we introduce PolyCLIP-T, a topology-guided multimodal framework that transforms variant selection from a classification problem into a geometric discovery task. By contrastively aligning DNA-sequence embeddings with functional annotations, PolyCLIP-T constructs a unified latent space in which the displacement between reference and alternate embeddings quantifies the molecular perturbation induced by each variant. Persistent homology then identifies stable topological components - coherent variant groups shared among affected relatives - that transcend single-variant scoring logic. Applied to six families with multi-morbid cancer, autoimmune and cardiovascular disease, PolyCLIP-T recovered non-coding and structural candidates overlooked by conventional pipelines and revealed pleiotropic networks spanning disease categories. This approach provides an interpretable, scalable solution for genome-first investigations of disorders driven by polygenic architectures that evade single-variant analysis. The framework was developed and benchmarked on deeply characterised familial cohorts selected for transgenerational multimorbidity; validation in larger, independent populations will be essential to establish its generalisability. An interactive web tool is freely available at https://www.polyclip-t.uma.es/.